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With the increased availability of high intensity time-of-flight neutron and synchrotron X-ray 
scattering sources that can access wide ranges of momentum transfer, the pair distribution function 
method has become a standard analysis technique for studying disorder of local coordination spheres 
and at intermediate atomic separations. In some cases, rational modeling of the total scattering 
data (Bragg and diffuse) becomes intractable with least-squares approaches and necessitates reverse 
Monte Carlo (RMC) simulations using large atomistic ensembles. However, the extraction of mean¬ 
ingful information from the resulting atomistic ensembles is challenging, especially at intermediate 
length scales. We use representational analysis to describe displacements of atoms in RMC ensem¬ 
bles from an ideal crystallographic structure in an approach analogous to tight-binding methods. 

Rewriting the displacements in terms of a local basis that is descriptive of the ideal crystallographic 
symmetry provides a robust approach to characterizing medium-range order (and disorder) and 
symmetry breaking in complex and disordered crystalline materials. This method enables the ex¬ 
traction of statistically relevant displacement modes (orientation, amplitude, and distribution) of 
the crystalline disorder and provides directly meaningful information in a locally symmetry-adapted 
basis set that is most descriptive of the crystal chemistry and physics. 


I. INTRODUCTION 

Achieving an atomistic description of solids continues 
to provide a challenge to the study of materials, especially 
as we learn that imperfections and disorder of crystals 
can give rise to the emergence of unexpected materials 
properties. For example, the multifunctional properties 
of the perovskite manganites can only be explained by 
understanding the relationships between the local and 
average structures mug. Therefore, we strive to fur¬ 
ther classify and quantify the nature of any local ordering 
(short-range order) that is patterned in a disordered fash¬ 
ion. Pair distribution function analysis of total scattering 
data has become a common technique for the characteri¬ 
zation of local distortions and disorder in crystals, as well 
as nanoparticle structures m- 

Modeling of atomistic structures - with emphasis on 
capturing the correct local structure - from experimen¬ 
tally derived atom-atom histograms poses a great chal¬ 
lenge, especially when the best description of the PDF 
has a short, finite correlation length (a domain) that be¬ 
comes averaged into a higher symmetry in the crystal¬ 
lographic structure. To obtain an atomistic description 
of such a model with these domains (where each domain 
consists of a few unit cells), simulations containing thou¬ 
sands of atoms can be used to model the total scattering 
data. By employing a large-scale simulation, the limita¬ 
tions from periodic boundary conditions are lifted, thus 
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allowing disordered aspects of the structure to either av¬ 
erage out into the Debye-Waller factor in the case of crys¬ 
talline disorder or lack any attributes of long-range order 
over the range of data provided in reciprocal space (af¬ 
ter convolution with the finite size of the simulation) in 
order to describe amorphous solids [7HH1- However, anal¬ 
ysis of these large-scale atomistic ensembles containing 
1000’s of atoms has been non-trivial, both in the chal¬ 
lenge of extracting information relative to the average 
crystallographic structure and also in providing statisti¬ 
cally meaningful information; there are typically many 
more free parameters in these simulations than indepen¬ 
dent observations (i.e., data). 

Herein, we develop a systematic approach for analyz¬ 
ing the disorder in large atomistic simulations of complex 
crystal structures using representational analysis. The 
determination of crystallographic superstructures result¬ 
ing from displacive distortions via symmetry-mode anal¬ 
ysis of a statistical distribution of ensembles has proven 
to be very powerful (cf., WO3 and LaMnOs) [ID] . An¬ 
other similar approach, but coupled to a different anal¬ 
ysis, has also made it possible to extract phonon disper¬ 
sions from powder diffraction data [mm. Here, we use 
a variation of this technique adapted to understand lo¬ 
cal structural variations by projecting displacements of 
atoms from their average crystallographic sites in atom¬ 
istic ensembles onto a tight-binding-like basis formed 
from the symmetry-adapted |J4] modes of a single unit 
cell, as depicted in Figure [l] - we define these modes 
as tight-binding modes. When displacements from an 
ideal crystallographic site are projected onto this locally 
symmetry-adapted basis, the disorder can be quantified 
and statistically analyzed to determine the frequency of 
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FIG. 1. Schematic illustration of the workflow for generating 
a large ensemble simulation using reverse Monte Carlo simu¬ 
lations, which can be folded back onto a small unit cell or the 
atomic displacements can be projected onto the tight-binding 
modes of small unit cell. 


specific displacement magnitudes and orientations. This 
manuscript outlines the analytical method and presents 
two illustrative applications of the method: the obser¬ 
vation of a trigonal distortion in BaTiOs at room tem¬ 
perature and the identification of the local displacement 
modes in the charge-ice pyrochlore, Bi 2 Ti 2 O 7 . More 
broadly, our approach is equally important for the anal¬ 
ysis of experimental diffraction and scattering data HU- 
Et ab initio and force-field-based simulations nans], 
and combinations of the two mm- Furthermore, this 
approach provides a common language and representa¬ 
tion for bridging experiment- and theory-derived models. 


II. METHOD 

A. Introduction to Total Scattering Methods 

The analytical method described here operates on an 
ensemble of atoms that can be described as a enlarged 
big-box generated from small crystallographic unit cells. 
The atom positions need not sit on precisely ordered lat¬ 
tice sites; however, upon back-folding the big-box ensem¬ 
ble onto the parent unit cell, the average atom positions 
should project close to particular lattice sites, each with a 
position distribution resembling something like a Debye- 
Waller factor (z.e., the model may be paracrystalline). 
This method is agnostic to how the models are gener¬ 
ated; the authors refer the reader to Refs, 
for a description of modeling total scattering data. 

Here, we use “total scattering” to refer to the scat¬ 
tering of X-rays or neutrons that describes the structure 
factor of the crystallographic symmetry (diffraction from 
periodically ordered components) and the diffuse scatter¬ 


ing that can arise from displacements of atoms from their 
ideal lattice points, including displacements from thermal 
motion and static disorder in the crystal [3]. If the total 
scattering structure factor, S(Q ), is measured to suffi¬ 
ciently high momentum transfer ( Qmax > 15 A -1 ), one 
can numerically take a sine Fourier transform to convert 
S(Q ) into the reduced pair distribution function (PDF), 
G(r), 

G(r) =4Trrp 0 [g(r) - 1] (1) 

OO 

Q[S(Q) - 1] sin (Qr)dQ, 

where po is the average number density of the material 
and g(r) is the atomic pair distribution function. The 
atomic pair distribution function, g(r ), is a direct mea¬ 
sure of the relative positions of atoms in a solid, that 
is, an experimentally-accessible real-space histogram of 
all atom-atom separations in the solid (of periodically 
ordered and disordered atoms). Because of the crystallo¬ 
graphic phase problem, without the use of isotopic label¬ 
ing or anomalous scattering, it is not possible to directly 
assign peaks in the PDF to specific atoms; therefore, 
atomic-scale modeling must be used to make assignments 
to individual peaks. 

“Small-box” models, which allow extraction of bond 
lengths and a description of the thermal motion (i.e., 
Debye-Waller factors), can be obtained from least- 
squares (LS) optimization of a crystallographic unit cell, 
or some small variant thereof to the experimental PDF 
using the software, PDFgui. %mmm LS optimization 
is susceptible to finding local minima in the goodness- 
of-fit and is numerically cumbersome when the model 
contains many degrees of freedom, as applicable here. 
Additionally, these short-range ordered models often fail 
to provide an accurate description of the crystallographic 
observations |26H28] . 

A complimentary approach to extract atomistic config¬ 
urations from the PDF is to simultaneously model both 
the crystallographic structure factor and the PDF by 
employing a “large-box” simulation of the total scatter¬ 
ing data. A Reverse Monte Carlo (RMC) algorithm can 
be used to find atomistic configurations of the ensemble 
consistent with both the experimentally determined G(r) 
and S(Q) [221 [23]. 

B. Coordinate Transform and Decomposition 

The goal of this method is to define the atomic config¬ 
urations of an ensemble as displacements from the ideal 
crystallographic positions. Here, we define a “large-box” 
as a M x x M y x M z enlargement of the crystallographic 
unit cell to form an atomistic ensemble, but no attempt 
is made to constrain the symmetry between atoms, either 
within the subcells or the “large-box”. The simplest such 
basis is to simply write down displacement vectors, in 
Cartesian or lattice coordinates, for each atom in the en¬ 
semble. Each atom within the crystallographic unit cell, 
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i, has a unique position defined by a vector, x^ n . The 
vector, R n , describes the spatial vector between each unit 
cell, n, within the ensemble. Each atom can be mapped 
as a displacement from its ideal position in the crystal¬ 
lographic unit cell, by u^ n = x^ n — where the 
values x^ n are often determined from a traditional crys¬ 
tallographic analysis (Rietveld analysis or single-crystal 
structural refinement). Such a representation is shown 
schematically for a simple two dimensional toy model, a 
2 x 1 “big-box” built from a crystallographic unit cell 
with two atoms and C 4 symmetry, in Figure [2|a). While 
straightforward to compute, this basis (the displacement 
vectors) lacks any connection to the symmetries that are 
present, locally or on average, and are thus difficult to 
interpret. A more refined approach is to rewrite the dis¬ 
placements in terms of the normal modes of the crys¬ 
tallographic structure, with amplitudes and phases for 
every mode at every wavevector in the Brillouin zone (as 
determined by the point symmetry of each wavevector). 
This normal-mode basis provides physical insight because 
the atomic displacements are mapped onto symmetry- 
defined motions away from their ideal positions, and cor¬ 
relations between unit cells are captured. There is, how¬ 
ever, an even better choice of basis that keeps many of 
the advantages of the classic normal mode approach but 
retains physical insight into the local symmetry changes. 

First, the local tight-binding (z.e., locally symmetry- 
adapted) modes are identified. This is accomplished by 
rewriting all possible atomic motions within a single unit 
cell into motions consistent with the point symmetry of 
the crystal at the Brillouin zone center, k = (0,0,0). 
Each motion (or mode) can be labelled according to the 
irreducible representation that it transforms under in the 
point symmetry group and is described by a set of basis 
vectors describing the actual atomic motions. Identifi¬ 
cation of these tight-binding modes is straightforward: 
basis vectors spanning each irreducible representations 
for each space group have been tabulated by Kovalev 
m, or can be computed by various crystallographic 
tools, including karep m, SARA h EQ , BASIREPS [32] , 
the Bilbao Crystallographic Server (Symmetry Adapted 
Modes) [33H36] . or the ISOTROPY software suite [37] . 
The input for these tools are the crystallographic space 
group and atom positions of the small (crystallographic 
average) unit cell as one would derive from Rietveld (or 
other suitable crystallographic) analysis. 

These tight-binding modes provide an orthonormal 
and local basis for describing all possible motional de¬ 
grees of freedom within a single unit cell, and are anal¬ 
ogous to normal vibrational modes of a molecular sys¬ 
tem. To retain this physical intuition but capture the 
degrees of freedom of a “large-box” atomistic ensemble, 
we adopt a technique analogous to tight-binding methods 
in electronic structure calculations [38] and write down 
modes at a non-zero wavevector in terms of these local 
basis functions that we define at the Brillouin zone center 
(a la atomic orbitals), with appropriate phase factors to 
describe correlations between crystallographic unit cells 
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FIG. 2. Schematic illustration of the coordinate transforma¬ 
tion from atomic displacements into tight-binding basis based 
on locally symmetry-adapted modes. 


in an ensemble. Specifically, we define the spatial cor¬ 
relations between unit cells within the ensemble with a 
quantized, reciprocal wavevector, k = 2 tt/R. The vector 
spans the indices, k x = (27rn x /M x ,27rn v /M v ,27rn z /M z ) 
for all n x = 0, 1 , • • • (M x - 1), n y = 0,1, • • • (M y - 1), and 
n z = 0,1, • • • (M z — 1); in other words, the wavevectors 
are in steps of 2tt/M along each direction. For mathe¬ 
matical convenience, we define all values for k as positive. 
The amplitude of a tight-binding mode, Tj^\(k), with the 
associated phase factor described by the reciprocal-space 
wavevector, &, is defined by: 

= Y y, Uj,n ■ ■ R n ), (2) 

n i 

where i runs over all atoms in the crystallographic unit 
cell, and n runs over all unit cells contained within the en¬ 
semble. The vector, i?, points to the n th crystallographic 
unit cell in the ensemble. The values, are the vec¬ 

torial contribution of atom i to the mode described by the 
(j, A) pair. The vectorial contributions can span multiple 
atoms, as pertaining to the crystallographic multiplicity 
of the particular site in the original crystallographic unit 
cell. The index, j, specifies each set of modes that to¬ 
gether transform as an irreducible representation of the 
point group; A is equal to the dimensionality of the cor¬ 
responding irreducible representation and runs over all 
modes in the set. Together, there are 3N distinct j, A 
pairs, or tightly-binding modes, where N is the number 
of atoms in the small, crystallographic cell. 

There is no index k on ^ J us f like there is no 
wavevector dependence to atomic orbitals in the clas- 
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sic tight-binding electronic structure approach, because 
all wavevector dependences are explicitly included in the 
phase factors. Further, note that to retain all degrees of 
freedom we allow the amplitudes of each tight-binding 
mode to be independent of all others, even if symmetry 
would constrain them (i.e., because one irreducible rep¬ 
resentation may be spanned by multiple modes). This 
allows for us to consider, but not enforce, symmetry in 
describing the “large-box” atomistic ensembles. Stated 
differently, the projection is only a change of basis; all 
3N — 6 degrees of freedom (for an ensemble of N atoms) 
are retained (omitting the 3 translational and 3 rotational 
degrees of freedom), and the exact atomistic ensemble 
can be reconstructed by the inverse of : 

-% = £££r A exp (ik -Rn)- (3) 

k 3 A 

This method, as applied to the toy model, is shown in 
Figure |2|b). We note that this is distinct from typical 
crystallographic order parameter analysis HMniEaEsi, 
in which the constraints of the parent crystallographic 
symmetry are preserved and the primary interoperable 
variables are the order parameter amplitudes, thus pro¬ 
viding one number for a pair of basis vectors that describe 
a displacement transforming as a two-dimensional irre¬ 
ducible representation, vs. two numbers in our approach. 
We retain all possible degrees of freedom. 


C. Continuous symmetry measures 

When using our tight-binding modes, we can determine 
the activity of the mode and deviation of the ensemble 
from the crystallographic symmetry not just from the 
mode amplitude, but also from its mean-squared devia¬ 
tion (MSD) from an ensemble operated on by a symme¬ 
try operation of the parent crystallographic space group. 
There are at least two distinct types of continuous sym¬ 
metry measures (as developed by D. Avnir and coworkers 
mm) that we characterize here. First, the global ac¬ 
tivity of a single tight-binding mode (j, consisting of one 
or more individual A modes depending on the dimension¬ 
ality of the corresponding irreducible representation) can 
be quantified as the mean-squared deviation (MSD) be¬ 
tween the (|r i>A (ife)|) amplitudes and the new amplitude 
coefficients, \TGj,x(ky\, following application of a sym¬ 
metry operation, G, of the crystallographic space group: 

s G,j = (£(£^4)? —Ei r ^(fc)l 2 ) 2 ) 1/2 - (4) 

k x x 

For purely symmetry-conserving displacements, the MSD 
should be zero. Here, it is critical to combine the squared 
amplitudes of all individual modes that together trans¬ 
form as a single multidimensional irreducible represen¬ 
tation (the innermost sums) because the amplitude of 


individual modes can be varied simply by changing the 
choice of basis vectors within that mode set. The sum 
over all wavevectors is justified to identify local symmetry 
changes because it corresponds to summing up contribu¬ 
tions derived from a single local tight-binding mode (a la 
atomic orbital), and is exact in the molecular limit. The 
final square root is provided for convenience to make the 
magnitude of sqj more physically interpretable. 

The related MSD, not broken down by individual mode 
sets, is similarly simple to calculate: 

s g = (££(£|r G j,A(fc) / l 2 -£l r i ) AWI 2 ) 2 ) 1/2 - 

j k x x 

(5) 

Where again the final square root is provided for conve¬ 
nience. 

The second type of deviations from the parent space 
group that can be identified are distortions that do not 
retain an equivalence of mode amplitudes within a sin¬ 
gle mode set that transforms as a multidimensional ir¬ 
reducible representation. To illustrate this, consider a 
single box with C 4 symmetry and an atom in the center 
displaced along the diagonal direction (Figure [ 3 ]). Pro¬ 
jected onto the two basis vectors T 1 and T 2 , which to¬ 
gether span the two-dimensional, E , irreducible represen¬ 
tation in the corresponding point group, the amplitudes 
along each basis are initially equal. As the Euler angle 
that defines the absolute orientation of the basis vectors is 
varied, the intensity of T 2 reduces while T 1 rises until T 1 
is collinear with the atom displacement; this oscillatory 
pattern continues the rest of the way. Note that the sum 
of the square amplitudes from the two contributions (Ti 
and T 2 ) is a constant (this is required as the magnitude of 
the displacement is not changing). However, across mul¬ 
tiple subcells or across multiple simulation runs, one can 
differentiate between random and ordered displacements. 
Let 0o be the initial angle of the displacement of the cen¬ 
ter atom. Different 0o’s correspond to phase shifts of the 
values of (and r|). If the 0o’s are completely random, 
then their average is a flat line as a function of Euler an¬ 
gle, with variances that are also flat [Figure j3^b)]. On 
the other hand, if the 0o’s are pinned to specific direc¬ 
tions, then only a subset of the phase shifts are present. 
This will often result in an average that is still flat as a 
function of Euler angle, e.g., if they are pinned every 90°, 
but the variances will no longer be uniform [Figure [3|c)]. 
This can be leveraged to determine whether the displace¬ 
ments are approximately random or fixed in some subset 
of orientations relative to the parent unit cell coordinate 
system. 

III. CASE STUDIES 

A. Trigonal displacements in tetragonal BaTiOs 

The ferroelectric ceramic BaTiOs at T = 298 K pro¬ 
vides an excellent example of a local distortion that av- 
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FIG. 3. Schematic illustration of the the rotation of the Eu¬ 
ler angle about the rotation axis of a particular symmetry 
element (C 4 ). (a) In a single box, the amplitudes of the dis¬ 
placements oscillate as a function of Euler angle, 0. (b) For 
random displacements across multiple subcells or boxes, the 
amplitude averages to a constant value as a function of 4> 
with a high and constant variance (denoted by error bars), 
(c) For displacements towards the four corners, the ampli¬ 
tude averages to a constant value as a function of 0; however, 
the variance oscillates as a function of <j) angle. 


erages out to a higher crystallographic symmetry in the 
unit cell. The average crystallographic symmetry deter¬ 
mined from Rietveld analysis is tetragonal, P4mm, which 
was used to define the tight-binding modes. However, the 
local bonding environment is significantly distorted and 
better described by the symmetry of the low-temperature 
R3m configuration [42H44] . as illustrated in Figure |4ja) 
and (b). While the R3m model provides a quantitative fit 
to the PDF within one unit cell, it does not provide quan¬ 
titative information of the medium-range order, such as 
information of the correlations between unit cells or the 
coherence length scale, even though such information can 
(and should) exist within the PDF. 

The experimental data used for this analysis were col¬ 
lected using the NPDF instrument (Lujan Neutron Scat¬ 



(a) Rietveld: PAmm 


(b) PDF: Ft 3m 



FIG. 4. (a) Room temperature PAmm crystal symmetry 

of BaTiOa with exaggerated displacements of the Ti and O 
positions to illustrate the ferroelectric dipole, (b) the pair dis¬ 
tribution function analysis reveals a local distortion present 
at room temperature that resembles the low-temperature 
R3m crystal structure [42H44] . illustrated here with exagger¬ 
ated Ti and O displacements, (c) Folded atomistic big-box 
ensemble generated from a RMC simulation, overlain with 
the anisotropic thermal ellipsoids determined from small-box 
modeling of the PDF (R3m structure). The (200), (020), and 
(002) planes are shown to illustrate the net displacement of 
the oxygen atom positions rather than the titanium and bar¬ 
ium atom positions. 


tering Center, Los Alamos National Laboratory) and 
were reanalyzed with adjusted relative absorption cor¬ 
rections (such that a scale factor was not needed to fit 
the intensity G(r)); the experimental details and original 
report of the experimental data can be found in Ref. 44 . 
The Bragg profile and PDF were used to constrain RMC 
simulations using the RMCprofile code, [22] as illustrated 
in Figure [lja) and (b). The simulation ensemble is a 
12x12x12 enlarged big-box ensemble of the tetragonal 
PAmm unit cell (8640 atoms) that was determined from 
Rietveld analysis. The ensembles were constrained by 
G(r ) (from 1 < r < 24 A) in addition to the Bragg profile 
from the 90° detector bank from NPDF (1.7 < Q < 15.7 
A -1 ; 3.7 > d > 0.4 A). In addition to hard-sphere cut¬ 
offs, a small penalty was applied to the simulations for 
breaking [TiOe] coordination in order to accelerate the 
simulations. 200 different simulations were performed 
from the same starting configuration in order to build 
statistics in the analysis. Each simulation ensemble can 
be back-folded into the unit cell; the atom positions fall 
within a cloud-like distribution centered around the av¬ 
erage crystallographic site [Figure Qc)]. 

Using the analysis method presented here, the atom 
positions were then decomposed into the tight-binding 
basis of the PAmm space group with a /c-mesh divided 
into 12 discrete steps along each x, y , and z direction 
with M x = M y — M z — 12. The irreducible repre¬ 
sentations and corresponding basis vectors for the tight- 
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FIG. 5. Visualization of selected tight-binding modes of 
BaTiOs in the PAmm space group. For Ti (16 Wyckoff po¬ 
sition), all three basis vectors are shown and demonstrate 
the retained 3 N degrees of freedom as the pair that together 
transform as E are allowed to have independent amplitudes. 
For the 02 site (2c Wyckoff position), the A\ and B\ modes 
each joins two oxygen positions, but 3 N degrees of freedom 
are retained for the two atoms generated from that Wyck¬ 
off position, noted by the six independent modes. While the 
E pairs will transform together if the local symmetry is also 
P4mm, all amplitudes are allowed to vary independently in 
this analysis. 


binding (locally symmetry-adapted) modes were identi¬ 
fied using the Bilbao Crystallographic Server (Symmetry 
Adapted Modes) [36] and are listed in Table [Tj some basis 
functions are graphically represented in Figure [5) 

For the analysis of a single ensemble of BaTiOs, the 
tight-binding mode amplitudes that describe displace¬ 
ments along the ferroelectric polarization are not very 
large (Ti A ± , 01 A u 02 A u 02 B u Table 0. This 
makes sense, since the average positions of the titanium 
and 02 atoms off-center along the elongated c axis di¬ 
rection (Table |) Hang. However, the displacements 
in the ab plane are significantly enlarged. This is graph¬ 
ically represented by the “point-cloud” distributions of 
the atom positions in Figure [4jc) that are overlain on 
top of the R3m unit cell used to describe the PDF in 
Ref. [44] . 

One problem with RMC simulations is that if the 
data are insufficiently resolving such that some atoms are 
poorly constrained, then the simulation atoms can wan¬ 
der away from their ideal positions. This would give the 
same graphical appearance as in Figure 0c); however, 
the quantitative data presented in Tablelll] shows that 


TABLE I. Basis vector components along each crystallo¬ 
graphic direction for BaTiOs, described in the PAmm space 
group setting, using the fractional atom coordinates, Ba: 
(0 0 0), Ti: (0.5 0.5 0.516), Ola: (0.5 0 0.487), 016: 
(0 0.5 0.487), 02: (0.5 0.5 0.978). 


Irrep 

Atom 

Basis vector components 
u\\a u\\b u\\c 

Ba Ai 


Ba 

0 

0 

1 

Ba E 

(i) 

Ba 

1 

0 

0 

Ba E 

(2) 

Ba 

0 

1 

0 

Ti Ai 


Ti 

0 

0 

1 

Ti E 

(1) 

Ti 

1 

0 

0 

Ti E 

(2) 

Ti 

0 

1 

0 

01 Ai 


Ola 

0 

0 

1 



016 

0 

0 

1 

01 Hi 


Ola 

0 

0 

-1 



016 

0 

0 

1 

01 E 

(1) 

Ola 

1 

0 

0 


(1) 

016 

0 

0 

0 

01 E 

(2) 

Ola 

0 

0 

0 


(2) 

016 

0 

1 

0 

01 E 

(3) 

Ola 

0 

0 

0 


(3) 

016 

1 

0 

0 

01 E 

(4) 

Ola 

0 

1 

0 


(4) 

016 

0 

0 

0 

02 Ai 


02 

0 

0 

1 

02 E 

(1) 

02 

1 

0 

0 

02 E 

(2) 

02 

0 

1 

0 


TABLE II. Amplitudes of the tight-binding modes (summed 
over all wave vectors) for the numerical average of 12 BaTiOs 
ensembles, a control simulation assuming PAmm local sym¬ 
metry, and a control simulation assuming R3m local symme¬ 
try. The values in parenthesis are the variance of the coeffi¬ 
cient. 



Tight-binding mode coefficients (A 2 ) 

Irrep 

Data 

PAmm control 

R3m control 

Ba A\ 

0.0034(1) 

0.0059(4) 

0.0026(3) 

Ba e\ 1) 

0.0127(4) 

0.0066(5) 

0.007(1) 

Ba E( 2) 

0.0130(3) 

0.0065(5) 

0.007(1) 

Ti A 1 

0.0089(3) 

0.0082(6) 

0.0024(3) 

Ti E( 1) 

0.0143(6) 

0.0077(6) 

0.009(1) 

Ti E(2) 

0.0128(5) 

0.0076(6) 

0.009(1) 

Ol A 1 

0.0044(1) 

0.0043(3) 

0.0031(4) 

Ol E( 1) 

0.0103(3) 

0.0068(5) 

0.008(1) 

Ol E( 2) 

0.0101(3) 

0.0068(5) 

0.008(1) 

02 A i 

0.0040(1) 

0.0045(3) 

0.0025(3) 

02 B x 

0.0039(1) 

0.0044(3) 

0.0024(3) 

02 E{ 1) 

0.0116(4) 

0.0076(6) 

0.009(1) 

02 E{2) 

0.0114(4) 

0.0076(6) 

0.009(1) 

02 E( 3) 

0.0090(3) 

0.0066(5) 

0.007(1) 

02 E( 4) 

0.0099(2) 

0.0066(5) 

0.007(1) 


these displacements are significant on average within an 
ensemble and that their variance is tightly defined even 
across 200 simulations. 

As a control, we performed RMC simulations con¬ 
strained by simulated pair distribution functions. In 
one case (PAmm control), we computed G(r ) from the 
PAmm crystal structure obtained by Rietveld analysis 
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FIG. 6. Euler angle analysis of a multidimensional tight- 
binding mode set that transforms as a multidimensional ir¬ 
reducible representation, (a) A random distribution of atom 
positions from the crystallographic location (dashed circle in 
cartoon) produces an equivalent variance of the mixing co¬ 
efficients (denoted with vertical bars) between the two basis 
vectors (ipi and ^ 2 ) when the basis vectors are rotated about 
the Euler angle, 0, parallel to the C 4 axis of the crystal struc¬ 
ture. (b) A clustering of positions at regular intervals, such 
as 7r/2, will produce the same mixing coefficients as in (a) 
for each tight-binding mode when averaged over all k and 
over all ensembles. However, a clustering of positions will 
yield a significant variance in the coefficients, denoted by the 
vertical bars, (c) The coefficients provided from the experi¬ 
mentally derived BaTiC >3 ensembles do not display significant 
differences when the basis vectors are rotated about the Euler 
angle. 


(convoluted with the appropriate instrumental resolution 
parameters, Qdamp and Qbroad ); the Bragg profile was 
the experimental Bragg profile. The simulated G(r) and 
Bragg profile were used to constrain 200 RMC simula¬ 
tions for analysis. For another control (R3m control), we 
took the reported R3m model determined from small-box 
modeling for the PDF (as reported in Ref. [44] ) and simu¬ 
lated G(r) from that structural model; the Bragg profile 
was the experimental Bragg profile). These then con¬ 
strained 90 independent RMC simulations for analysis. 
The PAmm control is a negative control that does not 
have additional displacements within the ab plane (be¬ 
yond thermal disorder modeled by a Debye-Waller fac¬ 
tor); the R3m control is a positive control for a known 
displacement in the ab plane coincident with thermal dis¬ 
order. The tight-binding mode coefficients resulting from 
analysis of the PAmm control simulations do not have a 
substantial anisotropy (Table |TT|); while there is a statis¬ 
tically significant increase in the coefficients of displace¬ 
ments in the ab plane, this may be biased from using 
the experimental Bragg peaks in conjunction with the 



0 (rad.) 


FIG. 7. Tight-binding mode coefficients as a function of Euler 
angle, 0, parallel to the C± axis of PAmm for (a) Ti, (b) 
0(1), (c), 0(2), and (d) Ba crystallographic sites, illustrating 
a constant variance as a function of <fi. 


simulated G(r). For the R3m control, there is a signif¬ 
icant and expected increase in displacements within the 
ab plane. This analysis informs us that the tight-binding 
mode amplitudes are capable of identifying aperiodic dis¬ 
placements when expected; however, the values of the co¬ 
efficients alone do not inform us to if particular symmetry 
operations are broken. 

With a local trigonal distortion, the R3ra-based model 
implies that there are specific vectors that the Ti dis¬ 
placements orient along - these are the vectors that point 
directly at the faces of the [TiOe] octahedra (i.e., the 
(111) directions, as referenced to the PAmm or Pm3m 
unit cells of BaTiOs). However, by looking at the graph¬ 
ical representation in Figure [4jc), it is impossible to 
tell if particular directions are preferred. Because the 
tight-binding modes with a set are mutually orthogonal 
and therefore yield locally orthogonal displacements, it is 
trivial to rotate the reference frame of the basis vectors 
and recompute their coefficient as a function of the Euler 
angle along the rotation axis of the multidimensional ir¬ 
reducible representation. In the PAmm description, this 
angle (</>) rotates around the 4-fold axis of the unit cell. 

In our analysis, we decomposed the atomic displace¬ 
ments into amplitudes of specific tight-binding modes as 
a function of rotation about the Euler angle, </> (Figure [6|. 
To illustrate this analysis, we employ two control simu¬ 
lations. Shown in Figure [6ja) is a simulation of displace¬ 
ments of Ti atoms around an approximately random dis¬ 
tribution of cf) angles. In Figure [6^b), we show a simu¬ 
lation with Ti atoms displaced at the same magnitude 
as in Figure |6|a), but the angles are constrained to be 
a random integer multiple of it/2 rad. Therefore, the Ti 
atoms cluster into four groups (akin to the (111) displace¬ 
ments). In both cases, the average coefficient of the tight- 
binding modes that together from a set and span a mul- 
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FIG. 8. Histograms of the mean-square displacements, 
summed over all k and mode sets (j) for 240 ensembles for 
each symmetry operation of PAmm , (a) the identity, E, (b) a 
vertical mirror along x, a v , ( c) a two-fold axis, 2 C 4 = C 2 , (d) 
a v + 2 C 4 , (e) a vertical mirror plane along xy , a v + Ca — cr xy , 
(f) the four-fold rotation axis, C 4 , and (g) 3Ch. 


tidimensional irreducible representation will not change 
as a function of </>, since the Ti atoms are displaced from 
the center by the same distance. However, the variance 
between tight-binding mode amplitudes [E( 1) vs. E( 2 )] 
will be distinct for each Euler angle (c/., Figure 3l). For 
the completely random distribution in Figure [6 v a), the 
coefficient multiplying T 1 of the E irreducible representa¬ 
tion will continuously vary between 0 and the maximum 
value, as the basis vector is orthogonal and collinear with 
the atomic displacement; the second basis vector (T 2 ) will 
also vary the same amount, but its amplitude will be 7 t/ 2 
out-of-phase with T 1 . Therefore, each basis vector will 
have the same variance with </>, denoted by the error bars 
in Figure [ 6 ] 

As in Figure J 6 ^b), if the atom displacements cluster 
into groups, then the variation of basis vector coefficients 
will not be constant with (j). When 0 is 0 rad., such that 
T 1 is oriented along the a axis, then its mixing coefficient 
will be y/2 times the average value; the coefficient of T 2 
will be identical. Therefore, the difference in coefficients 
is zero. However, when orients one of the basis vectors 
directly at the clustered displacements, one coefficient is 
maximal, and the other is zero; this produces a large vari¬ 
ation in the basis vector amplitudes. In the experimental 
simulations, there does not appear to be explicit cluster¬ 
ing of the titanium atoms along particular displacement 



(b) 


1 




FIG. 9. Histograms of the mean-square displacements of the 
(a, c, e) Ti A\ and (b, d, f) Ti E mode sets, as summed over 
all k for each unique symmetry orientation: (a,b) a x mirror 
planes, (c,d) a xy mirror planes, and (e,f) C 4 rotation axis. 


vectors [Figure [ 6 jc)]. Looking at the variation in coeffi¬ 
cients for all atoms in the unit cell, depicted by the error 
bars in Figure [7| there does not appear to be any cluster¬ 
ing of displacements as a function of Euler angle. While 
the two-dimensional E irreducible representations for 02 
appear to exhibit a trend with 0 , the change in average 
value of the coefficient reflects the definition of the ba¬ 
sis vectors; the variations of the coefficients, as indicated 
by the error bars, do not change with <fi. This result is 
consistent for RMC simulations run for different times 
(as disorder tends to be artificially maximized for longer 
simulation runs). 

While a variation of coefficients with Euler angle can 
indicate clustering of displacements described by a mul¬ 
tidimensional irreducible representation, it does not pro¬ 
vide any indication if the degeneracy-inducing symme¬ 
try operation is broken. To find broken degeneracies, we 
turn to continuous symmetry measures as defined in the 
methods. For BaTiOs, we compute the mean-square de¬ 
viation (MSD) for each generated symmetry operation 
of the crystallographic space group (PAmm). With four 
symmetry operations (E, cr v , C 2 , and C 4 ), there are a to¬ 
tal of 8 symmetry-related atoms that are generated from 
a general position; therefore, we test all unique combina¬ 
tions of these operations (each combination that gener¬ 
ates one of the general positions). 

The MSDs illustrate that the atomic displacements 
in the ensembles show the highest deviation away from 
four-fold rotation symmetry element. Histograms of all 
MSDs computed for BaTiOs (summed over all k and ir¬ 
reducible representations) are illustrated in Figure [ 8 ] for 
each symmetry operation. The histograms for related 
symmetry elements cluster together; those combinations 
that equate to a 4-fold rotation have the most significant 
MSD [Figure [ 8 ](g,h)] , followed by mirror planes parallel 
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to the {110} planes, then mirror planes parallel to the 
{100} planes. 

To probe which irreducible representations are most 
symmetry conserving, histograms of MSDs summed over 
all k for each irreducible representation are shown in 
Figure [9] the histograms bin together MSDs computed 
for the equivalent symmetry operations shown on the 
right. The histograms for the Ti A\ irreducible repre¬ 
sentation [Figure |9|a,c,e)] show tightly-grouped and low¬ 
valued MSDs indicating that the vertical Ti displace¬ 
ments tend to preserve the PAmm symmetry operations. 
However, the displacements that project onto the Ti E 
irreducible representation tends to break symmetry oper¬ 
ations, as previously inferred. The four-fold rotation axis 
appears to be the symmetry operation most broken, as 
naively expected from the small-box trigonal model illus¬ 
trated in Figure [4][b), which does have a vertical mirror 
plane parallel to the (110) plane. 

The analyses presented here for BaTiC >3 provide results 
that are sufficiently simple as to be easily compared with 
small-box models of BaTiC^. The use of RMC simula¬ 
tions allows one to extract a single, statistically relevant 
model of the atom positions that describes both the data 
regarding local atom separations (the PDF) and average 
crystallographic symmetry (Bragg profile). For BaTiOs, 
the coefficients of the tight-binding modes and their spa¬ 
tial dependence reveal the presence of a significant dis¬ 
tortion from the PAmm crystallographic symmetry. The 
resulting ensemble reveals that the atom positions are 
mostly displaced in the ab plane, which closely resem¬ 
bles the low-temperature R3m crystal structure. This 
example illustrates how such an analysis may be per¬ 
formed on materials with more complexity, both in terms 
of the crystal structure and the crystalline disorder, as 
described in the next section. 


B. Correlated Oxygen and Bismuth Displacements 
in Bi2Ti207 

The analysis methods presented here are generally ap¬ 
plicable to materials with more complex structures. The 
“charge-ice” pyrochlore oxide, Bi 2 Ti 2 07 , has a large unit 
cell that contains 88 atoms; direct inspection of en¬ 
sembles becomes prohibitive with this many degrees-of- 
freedom in a single crystallographic unit cell [37] . In 
Bi 2 Ti 2 07 , there is extensive disorder of the Bi sublattice, 
as attributed to stereochemical activity of the lone pair 
- derived from the [Xe]5d 10 6,s 2 electron configuration of 
Bi(III) - on a geometrically frustrated lattice. The geom¬ 
etry of the diamond lattice prevents long-range ordering 
of the dipoles, in a manner related to Pauling’s ice rules 
[48] . Previously, RMC simulations of total neutron scat¬ 
tering have been used to gain an atomistic representa¬ 
tion of the static Bi displacements, which form a toroidal 
distribution of Bi atoms that encircle the ideal crystal¬ 
lographic site. Furthermore, the O' atoms (8a Wyckoff 
site) are connected to the non-spherically distributed Bi 


BiE " 

«V 

Bi r -< 2 > S 3* 

FIG. 10. (a) Visualization of the tight-binding modes for the 
different Bi irreducible representations. The E u , T 2 u , and 
Ti u (l) mode sets have the largest amplitude, which describe 
the displacements that generate toroidal distribution of Bi 
positions and agree with the predicted displacement modes 
of E u and T\ u symmetry from ab initio density functional 
theory calculations. mi- 


atoms and therefore become pulled off of their ideal crys¬ 
tallographic site into tetrahedral volumes. The original 
report, experimental data, and experimental details can 
be found in Ref. [15 . The crystallographic Bi 2 Ti 2 07 unit 
cell is described by the Fd3m space group, which defines 
the irreducible representations and tight-binding modes 
used here. 

By rewriting the atomic displacements in terms of 
the tight-binding modes, a straightforward examination 
of their amplitudes reveals several characteristics that 
lead to many of the same conclusions as presented in 
Ref. [I5j; these coefficients are tabulated in Table [TTT| (the 
coefficients are averaged across modes related by face¬ 
centering, over all fc, and across 320 distinct ensembles). 
Of the Bi modes (depicted in Figure [To|), those spanning 
the E u and T 2w irreducible representations generate dis¬ 
placements that reproduce the toroidal distribution of Bi 
positions observed in Ref. m and have the most signifi¬ 
cant amplitude. The modes spanning the A 2u representa¬ 
tion is orthogonal to the C oo rotational axis of the torus 
and has a small amplitude. The modes spanning the 
Ti u (l) and T\ u {2) representations have intermediate ori¬ 
entations and amplitudes. The decomposition of atomic 
displacements into tight-binding modes reproduces the 
physically meaningful and intuitive results presented in 
Ref. [15]; here, the averaging across many wavevectors 
and simulations identifies the robustness of these dis¬ 
placements. 

Furthermore, identification of these high-amplitude 
modes allows one to create a “small-box” model for a 
symmetry constrained refinement. In Refs. m and ESI, 
imaginary phonon modes were discovered at the Brillioun 
Zone center; the symmetry of the polarization eigenvec¬ 
tors belong to the Tin and E u irreducible representations. 
The tight-binding modes spanning these irreducible rep- 
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Bi E u 

0.068(1) A 2 



a ^b Bi A^ g 

0.01023(2) A 2 


FIG. 11. (a) Visualization of tight-binding modes correspond¬ 
ing to the Bi T 2 u modes of Bi 2 Ti 207 . (b) The Bi T^ u mode 
set shows the Bi atom displacing from its nominal 180° O 7 - 
Bi-O 7 angle, orthogonal to the linear axis; this is reflected in 
the small magnitude of the Bi A 2n -derived modes. 



FIG. 12. (a) The O K\ g representation carries the largest 

magnitude of all degrees of freedom for the 48/ oxygen atoms, 
(b) Vectors depicting the O A\ g mode surrounding an linear 
[0 7 -Bi-0 7 ] group illustrates how the oxygens move out of the 
plane in which the modes of the Bi T<i u representation act (as 
illustrated in Figure [TT]). While all vectors are illustrated, not 
all of the atoms shown are spanned by a single irreducible rep¬ 
resentation (i.e., those related by face-centering are included). 


resentations have high-amplitude coefficients in the anal¬ 
ysis performed here. Distortion of the Fd3m lattice along 
these polarization modes yields a small unit cell of Cm 
symmetry that provides an excellent description of G(r) 
for r < 3.5 A m, The agreement of the high-amplitude 
tight-binding modes with the theory-predicted distortion 
modes and small-box refinement illustrates another util¬ 
ity of this approach for unknown systems. 

Additionally, the representational analysis performed 
here suggests that there are correlated Bi as well as 
O, Bi and O displacements that are not immediately 
observed from direct inspection of the atomic displace¬ 
ments. While a correlation between the Bi and O' dis¬ 
placements was previously made m, the tight-binding 
mode amplitudes show that there are large displacements 
of Bi and O. The highest modes corresponding to the 48/ 
oxygen atom correspond to the O A\ g irreducible repre¬ 
sentation, which can be described as a subtle elongation 
and twisting of the [TiOe] octahedra [Figure [l2|A)]. This 
large displacement is also mirrored in the anisotropic 
atomic displacement parameter of the 48/ oxygen po¬ 
sition obtained from Rietveld analysis. A possible origin 
of the high amplitude of this distortion is illustrated in 
Figure p^b): the 48/ O atom positions form a hexagon 
encircling the linear O'-Bi-O' linkages in an orthogonal 
orientation. With significant Bi displacements, as in¬ 
dicated by the large amplitude of the Bi T^ u spanning 
modes, the oxygen atoms displace from their ideal posi¬ 
tions around the hexagon in order to accommodate the 
shifted Bi atoms. 

By comparing different simulation runs, it is possible 
to gauge the uncertainty in how distorted or ideal the 
connectivity is in different parts of the lattice. For ex¬ 
ample, the mode amplitudes corresponding to the Ti-0 
sublattice are shown in Figure [l3j The histograms show 
that the displacement amplitudes have a narrow distribu¬ 
tion across all length scales and between many simulation 
runs. Furthermore, the displacements appear to be rea¬ 
sonably isotropic, as consistent with thermal disorder of 
a cubic lattice. 

In contrast, the distribution of Bi atom displacements 



mean-square deviation (A 2 ) 

FIG. 13. Histograms of the mean-square displacements of 
the (a) Ti and (b) O mode sets, as summed over all k , all 
symmetry operations, and 320 ensembles illustrate that the 
Ti-0 sublattice is not disrupted. Only the O A\ g mode shows 
any distribution of the MSD. 


are varied (Figure [l4| ). The Bi A^ u mode distribution 
is comparable to the Ti-0 sublattice. However, many of 
the Bi displacement modes corresponding to multidimen¬ 
sional irreducible representations have high amplitudes 
and broad distributions, indicative of substantial static 
disorder in directions orthogonal to the linear O'-Bi-O 7 
bond axis. This strongly suggests that those displace¬ 
ment modes locally break the Fd3m symmetry of the 
crystal structure. 

In this analysis, the multidimensional irreducible rep¬ 
resentations are broken into their individual components 
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TABLE III. Tight-binding mode amplitudes obtained from averaging 320 RMC simulations of B^T^CR. Atom positions are 
defined as Bi: (0,0,0), Ti: (0.5, 0.5, 0.5), 01: (0.43114, 0.125, 0.125), 02: (0.125, 0.125, 0.125); amplitudes are also averaged 
across atoms that are related through face-centering. 


Irrep Coef. (A 2 ) 

Irrep Coef. (A 2 ) 

Irrep Coef. (A 2 ) 

Irrep Coef. (A 2 ) 

Irrep Coef. (A 2 ) 

Irrep Coef. (A 2 ) 

Bi A 2u 0.010(1) 

Bi E u ( 1) 0.067(8) 

Bi E u ( 2) 0.067(7) 

BiT 2u (l) 0.068(8) 
BiT 2w (2) 0.067(8) 
BiT 2u (3) 0.067(8) 
BiT iu l(l) 0.030(4) 
BiTi u l(2) 0.030(4) 

Bi Tiit 1(3) 0.029(3) 

Bi T lu 2(l) 0.049(6) 

Bi T lu 2(2) 0.049(6) 

Bi T iu 2(3) 0.048(6) 

Ti A 2u 0.003(4) 

Ti E u ( 1) 0.008(1) 

Ti E u { 2) 0.008(1) 

Ti T 2u (1) 0.008(1) 

Ti T 2u {2) 0.008(1) 

Ti T 2u (3) 0.008(1) 

Ti T\ u l(l) 0.005(1) 
TiT ln l(2) 0.005(1) 
TiT lu l(3) 0.005(1) 

Ti T lu 2(l) 0.006(1) 

Ti T lu 2(2) 0.006(1) 

Ti T lu 2(3) 0.006(1) 

Ol A lg 0.021(2) 

Ol A 2u 0.006(1) 

Ol E u 1 0.006(1) 

Ol E u 2 0.006(1) 

Ol E g 0.006(1) 

Ol E g 0.006(1) 

Ol T 2 u 1(1) 0.007(1) 

Ol T 2 u 1(2) 0.007(1) 

Ol T 2 u 1(3) 0.008(1) 

Ol T 2u 2(1) 0.007(1) 

Ol T 2u 2(2) 0.007(1) 

Ol T 2u 2(3) 0.007(1) 

Ol T 2g 1(1) 0.007(1) 

Ol T 2g l{2) 0.007(1) 

Ol T 2ff l(3) 0.007(1) 

Ol T 2g 2(1) 0.006(1) 

Ol T 2g 2(2) 0.006(1) 

Ol T 2g 2(3) 0.006(1) 

Ol T 2g 3(1) 0.007(1) 

Ol T 29 3(2) 0.007(1) 

Ol T 2g 3(3) 0.007(1) 

Ol T iu l(l) 0.008(1) 

Ol T lu l(2) 0.008(1) 

Ol T iu l(3) 0.008(1) 

Ol T lu 2(l) 0.006(1) 

Ol T lu 2(2) 0.006(1) 

Ol T iu 2(3) 0.006(1) 

Ol T lu 3(l) 0.007(1) 

Ol T lu 3(2) 0.007(1) 

Ol T lu 3(3) 0.007(1) 

Ol T\ g l{l) 0.007(1) 

Ol Ti g l(2) 0.007(1) 

Ol Ti s 1(3) 0.007(1) 

Ol T lfl 2(l) 0.007(1) 

Ol T ig 2(2) 0.007(1) 

Ol T ig 2(3) 0.007(1) 

02 T 2ff (l) 0.029 (3) 

02 T 2g (2) 0.029 (3) 

02 T 2g (3) 0.029 (3) 

02 T lu (l) 0.030 (3) 

02 T lu (2) 0.030 (3) 

02 T lu (3) 0.030 (4) 



mean-square deviation (A 2 ) 

FIG. 14. Histograms of the mean-square displacements of the 
Bi mode sets, as summed over all k , all symmetry operations, 
and 320 ensembles, illustrate that the E u , T^ u and Ti u (2) rep¬ 
resentations break the symmetry operations of FdSm ; those 
modes correspond to toroidal displacements shown in Fig¬ 
ure [10] 


(as to retain the total number of degrees of freedom); 
however, the values of each tight-binding mode are iden¬ 
tical in the case of B^T^CV To then identify if and 
by how much atomic displacements break the symmetry 
elements linking together tight-binding modes spanning 
a multidimensional irreducible representation, the con¬ 
tinuous symmetry measure of each irreducible represen¬ 
tation can be calculated. Figure p~5] contains histograms 
of the MSDs for each irreducible representation after op¬ 
eration on the simulation box by a specific symmetry 
operation (Eqn. [4|. The modes spanning the A^ u rep¬ 


resentation do not show any dependence on the symme¬ 
try operation; however, the modes corresponding to E u 
and T 2u representations do show a dependence on the 
operations. Specifically, the face-centering [+(^,0, |), 
+ (0, !)] and inversion (i) symmetry operations show 

the highest MSDs, as well as the broadest distributions 
suggesting that those symmetries are deviated by the 
largest magnitude and in the most ways. In future work, 
it will be informative to analyze the compatibility rela¬ 
tionships as the degeneracy of different modes changes as 

M (o,o,o). 

The crystal structure of B^T^Oy presents a very com¬ 
plex problem, as the unit cell contains 88 atoms, resulting 
in 264 degrees of freedom, or 264 distinct tight-binding 
modes, to describe all atom displacements. When trying 
to analyze a large ensemble simulation of this structure, 
analysis in Cartesian coordinates becomes unwieldy. De¬ 
composition of the structure into the crystallographically 
relevant local basis allows one to determine the highest 
amplitude disorder in the lattice, the distribution of am¬ 
plitudes, the direction of the atomic displacements caus¬ 
ing the disorder, how the disorder breaks specific sym¬ 
metry elements of the crystallographic space group, and 
by how much. 

IV. CONCLUSIONS 

The representational analysis of large atomistic ensem¬ 
bles generated from simulations (such as from reverse 
Monte Carlo simulations of total scattering data) us¬ 
ing a tight-binding basis derived from locally symmetry- 
adapted modes is a robust method that allows one to 
quantify disorder in the lattice. In many RMC simula¬ 
tions, the goal is often to characterize subtle deviations 
from the lattice - these types of displacements are subtle 
perturbations from a lattice that posses a modicum of 
moderately isotropic thermal disorder. Therefore, isola¬ 
tion and quantification of the disorder (i.e., of infrequent 
events) requires statistical analysis. By representing the 
disorder with respect to a local basis of the background 
signal (ie., symmetry-adapted modes of the crystallo- 
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FIG. 15. Histograms of the mean-square displacements of the Bi mode sets for each symmetry operation of FdSm , as summed 
over all k and 320 ensembles. The amplitude and breadth of the MSD’s show that the face centering [+(|,0, |), +(0, |)] 

and inversion (z) symmetry operations are not well followed by the E u , T^ u and Ti n (2) representations. 


graphic space group), displacements appear as a posi¬ 
tive signal, are amplified and can be statistically ana¬ 
lyzed. Additionally, the approach presented here permits 
a framework for analyzing other types of degrees of free¬ 
dom, such as occupational/compositional disorder (e.g., 
solid-solutions) or magnetism. Such a rigorous group the¬ 
oretical treatment is currently implemented in ISODIS- 
PLACE [39]. 
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